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Motivated by the energetic advantage of achieving coherent enhancement of effective spin- 
dependent interactions through approximate nesting, we propose specific forms of spin ordering, 
whose form varies over the Fermi surface, for the cuprate superconductors. Competing "spin-orbit" 
orderings involving order parameters in spatial d x 2_ y 2 and d xy waves at commensurate and in- 
commensurate wavevectors, phase separated in momentum space, support behavior suggestive of 
observed phenomena. The d xy spin-orbit fluctuation induces an effective interaction that favors 
d 2 ,2_j / 2-wave pairing, as required for the observed superconductivity. Anisotropic spin susceptibility 
is a crucial prediction of our mechanism. 



High T c superconductivity in cuprates has been an out- 
standing problem in condensed matter physics since it 
was first discovered in 1986. In the meantime, though no 
consensus theoretical understanding has developed, sev- 
eral striking features of the phenomenology have emerged 
clearly. These include the strongly 2-dimensional char- 
acter of the essentially new physics; its proximity to an- 
tiferromagnetic spin ordering; the anomalous character 
of the normal state, especially in the underdoped region, 
suggestive of emergent energy gaps; the acute sensitivity 
of the low-temperature state to doping and impurities; 
the d-wave character of the superconductivity; and the 
apparent uniqueness or near-uniqueness of cuprate layers 
in supporting this overall phenomenological profile. On 
the face of it, several of these features suggest that the 
origin of these phenomena will involve forms of ordering 
that involve spin and depend sensitively on the precise 
form of the lattice and the Fermi surface. 

There is a simple heuristic that appears to be broadly 
consistent with these indications. As is familiar from 
the BCS theory of superconductivity, the effect of weak 
attractive interactions can be amplified, and can lead 
to drastic qualitative effects, if there are many low- 
energy pairs sharing the same quantum numbers. Corre- 
lations among these pairs can then be arranged so that 
their interactions contribute coherently to lowering the 
energy. The BCS mechanism of superconductivity in- 
volves particle-particle and hole-hole correlations. In 
this context the existence of low-energy pairs with to- 
tal momentum zero, arising from time-reversed states 
with momenta (k, — k) both near the Fermi surface, is 
generic. By contrast spin-density-wave ordering, at the 
level of electron creation and destruction operators, in- 
volves particle-hole correlations. In this context one finds 
that many low-energy pairs sharing a common (lattice) 
momentum only for specially shaped ("nested") Fermi 
surfaces. Since the shape of the Fermi surface changes 
with doping, one might anticipate that at best a nest- 
ing condition would be approximately fulfilled at a spe- 
cific doping level. Our point of departure is to consider 
that the Fermi surface is not an end in itself, but a step 



toward constructing the ground state. If changing the 
pattern of occupied levels — effectively, engineering the 
Fermi surface — can encourage favorable coherence fac- 
tors, it might be favorable to make nesting persist. Real- 
izations of this possibility and the properties of the emer- 
gent states will, on the face of it, depend sensitively on 
details of the interactions, lattice structure, and doping 
level. Two dimensional antiferromagnets on a square (or 
nearly square) lattice near half filling, as in cuprate lay- 
ers, provide an especially favorable area for these ideas. 

Proposed Ordering At half filling antiferromagnetic 
(AF) spin ordering is observed. Electrons in the real 
materials may well be best described as strongly coupled 
and the spins as localized, but we shall construct our 
states heuristically by extrapolation from weak or inter- 
mediate coupling, anticipating that universal properties, 
specifically symmetry breaking patterns, might be suc- 
cessfully inferred. (Also, photoemission experiments 
can be interpreted as revealing a Fermi surface, even at 
quite small doping.) In that spirit AF ordering can be 
regarded as follows. On-site Coulomb repulsion induces, 
in the crossed channel, an attractive interaction between 
electrons and holes of opposite spin at momentum trans- 
fer 7r = (7r,7r) (modulo reciprocal lattice vectors). This 
makes it favorable to deform the effective Fermi surface 
into a diamond shape, which for half filling nests at tv 
(see Fig.^i), and allow the electrons to form spin-triplet 
particle-hole pairs, according to (c k+7r <7Ck) 7^ 0. Here, 

(7 = (<n, (T2, 0-3) are Pauli matrices, and Ck and are 
electron annihilation and creation operators, with two 
spin components united into a column: = (c^, c^). 

At finite doping no ansatz seems so uniquely com- 
pelling, but the possibility illustrated in Fig. ^5 is sugges- 
tive. The free Fermi surface has been deformed in two 
distinct ways, one operating near to and the other far 
from the zone diagonals. This geometry supports nest- 
ing for orderings with wavevector tv in the off-diagonal re- 
gion and with wavevectors slightly off tv, namely, ±Q/, = 
((1 =F S)tt,tt) and ±Q„ = (n, (1 =p S)ir), along the diago- 
nals. This proposal embodies a new phenomenon, phase 
separation in momentum space, that might find wider 
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FIG. 1: Phase separation in momentum space, (a) Incommen- 
surate ordering at S — 2 arcsin(|/i|/2t) Q- B is suggested by 
approximate nesting for a Fermi surface with representative 
values of nearest and next-nearest neighbor hopping, t and t' , 
respectively, (b) A re-arrangement of occupied states allows 
incommensurate nesting to occur together with tt nesting for 
finite doping. 



application. 

Order parameters with near-uniform s-wave structure 
do not easily accommodate such phase separation. It is 
more natural when zeroes of one condensate correspond 
to maxima of the other. This leads us to propose spin- 
orbit (SO) orders of the form: 0] 



lE(r 



(4«v)(4,c r ) + (4acr>) ■ (cl,ac r ) (4) 



plus an unimportant term proportional to the density 
operator. In this form the anticipated electron-hole at- 
traction is manifest. (At this level Hy, with V > 
0, also favors charge-orbit (CO) ordering, known as 
orbital-antiferromagnetism, staggered flux phase [5j, or 
DDW 6] ; a non-static version was proposed in 0j. We 
shall not discuss it further here.) Similarly, the next- 
nearest-neighbor Coulomb interaction favors both d xy SO 
and CO orders. 




<4 + ,^ k ) = iN A T£ 
(k € off-diagonal) 



(j = x, y; k € diagonal) 



B 



(1) 

(2) 



Here Na,b are constant, real vectors in spin space, 
and r k ' are orbital wave basis functions of lattice 



cos k y and 



group D^h, defined by r k = cosfc^ 
sin k x sin k y . (The lattice constant is set to a unit.) The 
d x 2_ y 2 order is purely imaginary, as required by hemitic- 
ity. Roughly speaking, it describes a state with micro- 
scopic spin currents flowing around each plaquette in real 
space. 

The nature of these spin-orbit orderings may be more 
transparent in real space: 
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iN A e v "' r T?_ t 
x [ C0S (£^)+(^y)]rf_ r , (3) 



where Yf„ = §[(£. 



y ^r" ,x — y $r" 



,x ~t~ $r" ,—x) (*£ 



y)} and r* = 
-y] . Note that 
these vanish for r = r'. 

At mean field level such spin-orbit orderings with 
nonzero Na and Nb are favored by the nearest and next 
nearest neighbor Coulomb repulsions, respectively. In- 
deed the nearest-neighbor Coulomb interaction Hv = 



FIG. 2: Schematic phase diagram for cuprate superconduc- 
tors. The doping dependence of the transition temperatures 
Tg Q ~ y and T^q results from the competition of the two SO 
orders in phase space (see Fig. 0d). The doping value xi 
for the onset of dSC is estimated as the level of doping above 
which the Fermi surface near (±7T, 0) and (0, ±7r) starts to ap- 
pear. To the right and from above, the dSC state is bounded 
by the d^-SO order, upon which it depends. 



Phase diagram We shall focus on the following com- 
peting orders, that we believe play major roles: AF, 
d x 2_ y 2-SO and d xy -SO, and d x 2_ y 2-w&ve superconduct- 
ing state (dSC). The two SO orders compete for particle- 
hole pairing states, plausibly with different outcomes in 
different domains, whose size depends on overall doping 
level. This can be understood by reference to Fig. 
With increasing density x of doped holes, the total area 
of the four triangles along the zone diagonals increases 
at first while the size of off-diagonal regions shrinks. For 
sufficient large doping, e.g., x±, the triangles are no long 
sustainable. This explains the trends of phase transition 
lines of Tg ~ y (for d x 2_ y 2-SO) and TgQ (for d xy -SO) as 
functions of doping. 

For low doping up to Xz, there are several competing 
spin orders (including AF) that are connected by first- 
order phase transitions. On general grounds, one expects 
that phase separation (in real space) may occur, plausibly 
in the form of stripes H,EJ). 
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Effective Theory We now propose an effective La- 
grangian for the unconventional normal state, based on 
the hypothesis of d x 2_ y 2-SO ordering: 

L = cl a (id t + n)c ra -H + V A y^ y (j) r -<f r (5) 



where <p r = i £ r , T^_ r , \c\ac r i — h.c] is a composite op- 
erator, Va > is assumed, and H collects the hopping 
terms. The interaction term contains part of the original 
next-nearest neighbor Coulomb interaction. 

A d x 2_ y 2-SO ordered state is characterized by a non- 
vanishing order parameter field </>, corresponding to the 
imaginary part of c\ac v i , pointing to (say) the 3-direction 
in spin space, 



) = o, 



(6) 



This correlation spontaneously breaks SU(2) spin sym- 
metry down to a U(l) corresponding to spin rotation 
about the 3-axis. Two Nambu-Goldstone bosons appear 
as gapless collective spin excitations, corresponding to 
transverse fluctuations of the order parameter field <j) r . 
We will call them orbital magnons. 

We can describe the low energy interactions of these 
collective modes using an effective Lagrangian. Follow- 
ing the standard technique [Io|. we isolate the Nambu- 
Goldstone part (two transverse spin fluctuations) in the 
electron field. Expressing the electron field as a local 
SU(2) spin rotation acting on a new fermion tp: 

c r = 1 + i ^S^ Sl ^ = Ur ^ , ( a , /3= i j2 ; no 3). (7) 

i + Cl 

This defines a new fermion ip ar that carries the full charge 
and the spin (T3 quantum numbers but does not carry 
transverse spin. C, 1 ' 2 parameterize the slowly-varying or- 
bital magnons. They are given by the 3-component <f> 
via: U$a a U r = Rfa b and = Rfdfo) + g r ) with 
a, b = 1,2,3, where g represents the (gapped) longitudi- 
nal spin-orbit fluctuation. With the above transforma- 
tions, we now have a prescription of deriving an effective 
Lagrangian for the orbital magnons and new fermions: 
L[V>, V^i C Q ! 9} from the Lagrangian 

The free part of the effective theory for the ^-fermion 



L i> = Ek { ^llidt - e(k) + /j]ipk 

+*|A k |[^ k+7r ^k - h.c.]} 



(8) 



where e k = —2t(cosk x + cosk y ) — At' cos k x cos k y and 
A k = A^cosfca, — cosfcy) with = -VaNa- Due to 
unit cell doubling, the fermion energy spectrum is split 
into two bands 



each being spin up and down degenerate. For weak 
d x 2_ y 2-SO order, the Fermi surface of the E~ band is 
near the diagonals of the Brillouin zone while the Fermi 
surface of the E + band is near the off-diagonal region 
(see Fig. |3J). Note that the fermion field ip, appropriate 
to describing the low-energy excitations, does not carry 
the full spin quantum numbers of the electron. This is a 
form of partial spin-charge separation. 




FIG. 3: The Fermi surface in d„2_ y 2-S0 state. J5 k = is 



shown as the dotted ellipses, inside which all states are empty. 
E£ = is shown as the boundaries of the four patches (in 
dark) which are completely occupied. 



The effective theory for d xy -SO can be constructed in 
the same manner from a model Lagrangian similar to J5J , 
with V A -> V B and $ r -> $ r = 2 £ r , rf_ r , [c\ac r , + h.c] 
(d xy -SO order parameter). The ordering wavevector is 
changed from 7r to Qh,v 

Origin of d-wave Superconductivity In the effective 
theory of d xy -SO, fermions are coupled to the longitu- 
dinal spin-orbit fluctuation g, 

Hint = ~2V B £ kq [lt + rf+J^aV-k-qJfti ■ (10) 

Broken spin symmetry implies that the g propagator has 
the form (ui 2 — S e {(\} 2 + ij(q, w) 2 ) -1 with £ g (q) having 
minima at q — > ±Qh.v and the damping rate 7(q, u>) — > 
as w — > 0. Given the coupling l(Tfj|) . we have calculated 
the effective interaction between fermions induced by the 
longitudinal SO fluctuation (Fig. HJ. In the low energy 
(static) limit, the spin-singlet component of the Cooper 
channel interaction is described by in the Hamiltonian 



Ht° = E kp M^Jp^iJpc "» ^ 



(11) 



5(e k + e k +7r) ± \ /j(ek 



(9) 



where v kp = 2g g V%\Tg + T^] 2 S~ 2 {k - p) > and the 
spin indices p,a =T, J,. Here g e is a positive normaliza- 
tion constant of the SO field. The interaction is repulsive 
in both s-wave and d xy -w&ve pairing. Note that t> k p, as 
function of k — p, peaks at the incommensurate wavevec- 
tors ±Q/i jt) (Fig. UJ. By examining a typical supercon- 
ducting gap equation, A k c = - £ p w kp A p c [(e p - p,) 2 + 
Ap 2 ] -1 / 2 , one finds nontrivial solution if A^ and A_ c 
have opposite sign. By the arguments similar to Ref . , 
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FIG. 4: Effective attraction between spin-broken electrons 
induced by longitudinal spin-orbit fluctuation g. The arrows 
indicate spin. The free electron fermi surface is plotted for 
simplicity. The off-diagonal regions for dSC or d x 2_ y 2-S0 
need not to change qualitatively from the state of d xy -SO. 



cal spins does not exhibit conventional long range or- 
der. A straightforward exercise confirms (S r ■ S r >) — 
-| l(4 c r')| 2 + |K4«fcr')| 2 -» as |r - r'| -> c». 

We have calculated the uniform, static spin susceptibil- 
ity for the d x 2_ y 2-SO state, using the effective Lagrangian 
|jHJ. We find that the susceptibility is anisotropic, 

x ab = -Mi/^iran+ra:)] (13) 

x[f + 2(A£A^-^|A k | 2 )/(i : ;+- J B k )] 

where f(E) = el iE +l , f'(E) — df/dE, and the ordering 
direction is assumed arbitrary (cf. Eq. |JSJ). Fig. [5] shows 
the anisotropy of susceptibility predicted by the mean 
field theory (without corrections due to orbital magnon 
scatterings). Measurement of x" and x ± could therefore 
provide important tests of our proposals. The result is 
quite different from the AF state, for which x" < x^- 



we conclude that the favored superconducting pairing is 
(4.2 -wave. Fermions near (±7r, 0) and (0, ±7r) will pair 
as: (ipkpip-ka) ~ ^perTk 7^ . The V'-fermion pair conden- 
sate is equivalent to an electron pair condensate 



j)(7 



(iprplpr'v) = (CrU*i<J 2 Ul,C r i) = € pa (c r p C r > a ) (12) 



for uniform Q T . Thus an exotic normal state can be asso- 
ciated with conventional electron pairing. 

The state of d x 2_ y 2-SO is different. A similar coupling 
involving the longitudinal SO fluctuation in the d x 2_ y 2 
leads to an effective interaction that disfavors the d x i_ y i- 
SC pairing. 

d xy -SO order does not forbid the existence of Fermi 
nodal points in the dSC state. It splits the quasiparticle 
energy band around the four triangles along the zone 
diagonals (Fig. ^) , with Fermi surfaces being carved out 
from the lowest possible band. 

Incommensurate spin excitations The translation^ 
symmetry broken by the SO orderings leads to multi- 
bands of orbital magnons in a reduced zone scheme, each 
having the degeneracy of two transverse spin modes. The 
orbital magnons associated with the d xy -SO are gapless 
at the four incommensurate momenta iQh and ±Qt,. 
With coexisting d xy and d x 2_ y 2 SO orders, there will be 
another band that is gapless at 7r. Dynamical spin-spin 
correlation functions are affected by particle-hole pair 
mixing with two orbital magnon channels. Combining 
magnons from d xy and d x 2_ y 2 can yield two branches of 
spin excitation resonance at finite energies, one dispers- 
ing upward and another downward from momentum tt. 
The two branches cross at 7r, with zero gap energy. This 
structure seems qualitatively consistent with the inelastic 
neutron scatte ring experiments in YBCO (e.g, as in D. 
Reznik, et al. jl2| and references therein; see also [l3j|L 
A full investigation is in progress. 

Anisotropic spin susceptibility The state of SO or- 
dering has (S T ) = 0, so the static correlation of lo- 
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FIG. 5: Anisotropic spin susceptibility in the d x 2_ y 2-SO state 
at different temperatures. Energy units are chosen such that 
t = 1. Other parameters: t' = —0.45 and fi = —1.2 (chem- 
ical potential). Integration was performed numerically using 
a mesh of 500 x 500 in each quarter of the Brillouin zone. 
x"' x denote the susceptibilities parallel or perpendicular to 
the direction of the d T 2_„2 SO. 
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